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Accurate  prediction  of  variations  in  the  height  of 
atmospheric  inversions  are  required  for  optimum  utilization 
of  modern  weaponry.  A  marine  atmospheric  boundary  layer 
model  was  used  to  explore  the  influence  of  land  on 
predictions  of  the  inversion  heights.  Data  for  this  evalua¬ 
tion  were  obtained  from  rawinsonde  observations  taken  at 
Vandenberg  APB,  California  for  each  season.  Vertical  veloc¬ 
ities  fcr  these  data  were  obtained  by  vertical  integration 
of  the  moisture  budget  equation  (Q-method)  .  Results  of  the 
model  output  were  compared  graphically  and  by  use  of  RMS 
error  statistics.  These  results  were  very  poor  and  revealed 
that  strong  land  influences  on  diurnal  variations  yield 
inversion  height  fluctuations  which  the  model  cannot  repro¬ 
duce.  further  development  of  the  model  is  necessary  to 
overcome  problems  caused  by  the  land  influence. 
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I.  INljjOEgCTION 


Historically,  advanced  detection  and  surprise  attacks 
have  teen  an  area  cf  interest  to  the  military  community. 
Knowledge  about  the  effects  cf  atmospheric  boundary  layer 
(A3L)  inversions  on  electromagnetic  signals  is  well  known 
and  documented.  Advance  knowledge  of  inversion  changes  can 
have  a  direct  effect  on  military  planning  conducted  by  ail 
services. 

Kith  technolcgical  advances,  our  weaponry  has  fceccme 
highly  sophisticated  (e.g.,  the  development  of  "smart” 
tombs)  and  dependant  on  electromagnetic  propagation  for 
guidance  ccntrol.  The  guidance  components  for  these  weapons 
depend  cn  many  different  wavelengths,  ranging  from  the 
visible  tc  the  microwave,  dependent  on  the  specific  use  of 
the  weapon.  The  weapons  depend  on  the  transmission  and 
reception  cf  these  signals  to  locate  and  follow  targets. 
Cifficulties  arise  when  atmospheric  conditions  are  such  that 
signal  performance  is  distorted  by  refraction.  Only  these 
areas  which  affect  propagation  of  radar  signals  will  be 
discussed  in  this  thesis.  Other  propagation  types,  such  as 
optical,  and  the  atmospheric  parameters  affecting  th®m  will 
rot  be  addressed. 

Signal  propagation  is  dependent  on  the  index  of  refrac¬ 
tion  which  is  a  function  of  temperature  and  humidity. 
Eefracticn  cf  a  signal  occurs  as  that  signal  moves  through 
vertical  gradients  cf  temperature  and  humidity.  Atmospheric 
conditions "which  enhance  refraction  are  inversions,  which 
cause  strong  gradients  of  temperature  and  humidity,  or  rapid 
changes  cf  humidity  near  the  surface.  Occurrence,  height 
and  intensity  of  inversions  are  strongly  dependent  on  local 
and  synoptic  atmospheric  conditions.  The  importance  of 


predicting  propagation  is  obvious  but  to  do  sc  with  accuracy 
first  requires  accurate  prediction  of  the  met ecrcicgrcai 
parameters  affecting  it.  Of  obvious  import  a  r.ce  is  the 
prediction  of  changes  in  the  temperature  and  humidity 
profiles,  both  of  which  are  highly  variable  and  can 
completely  change  their  vertical  profile  within  a  matter  of 
hours.  Other  parameters  such  as  vertical  velocity  and 
surface  wind  speeds  will  be  discussed  later. 

The  ability  to  forecast  inversions  and  how  they  may  vary 
directly  influences  the  effectiveness  of  modern  weaponry. 
It  is  readily  apparent  that  improvements  in  this  forecasting 
ability  is  cf  paramount  importance  in  making  military  plans 
which  involve  choosing  the  type  weapon  best  suited  for  a 
particular  mission.  As  long  as  there  are  functions  depen¬ 
dent  cn  signal  transmissions  there  will  remain  a  need  to 
increase  the  ability  tc  forecast  the  atmospheric  parameters 
which  contribute  to  propagation. 

The  model  used  tc  investigate  inversion  heights  is  a 
marine  atmospheric  boundary  layer  (HABL)  model  which  has 
been  modified  for  use  at  tb%  Naval  Postgraduate  School 
(NPS )  .  The  factors  incorporated  into  this  model  include 
potential  temperature  and  specific  humidity  within  the 
well-mixed  region  and  their  gradients  at  the  base  of  the 
inversion.  Use  of  thesa  parameters  make  this  model  appli¬ 
cable  to  studying  inversion  changes.  Actual  discussion  cf 
the  medal  will  be  provided  in  Chapter  II. 

Another  area  where  advance  knowledge  of  ABL  inversion 
changes  vculd  be  useful  is  in  single-station  forecasting. 
The  effects  of  chances  in  the  inversion  are  found  not  only 
to  effect  electronic  transmissions  but  also  to  indicate 
possible  changes  to  current  weather  condirions.  Therefore, 
any  improvements  in  cur  ability  to  predict  these  changes 
will  be  important  tc  land-based  and  sea-based  operations, 
since  improvement  in  this  area  cf  prediction  also  will  be 
reflected  ir.  the  daily  weather  forecast. 
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This  thesis  completes  an  analysis  on  four  time  periods 
cf  data  obtained  from  Vandenberg  AFB,  California.  These 
data  periods  each  cover  a  different  season  and  vary  fro*  two 
to  eight  days  in  length.  In  conducting  this  examination 
vertical  velocities  fcr  Vandentarg  AFB  were  calculated  from 
their  rawinsonde  data.  Each  days  data  were  run  through  the 
NPS  S1EL  prediction  model  tc  establish  the  forecast  trend 
for  that  period.  Eesults  from  the  model’s  forecasts  were 
compared  to  data  from  later  reports  for  correlation  ar.d  to 
determine  accuracy  cf  the  model.  Successive  runs  of  the 
model  were  made  until  each  data  set  has  been  completed  and 
determination  as  to  how  well  the  model  predicts  inversion 
heights  in  each  24  hcur  forecast. 

One  cf  the  time  periods  of  data  chosen  for  study  is,  by 
choice,  coincident  with  a  data  pariod  used  by  Gleason 
(1982).  His  study  was  directed  at  effectiveness  of  the 
model  fcr  mere  of  a  pure  marine  location  and  for  this  reason 
data  from  San  Nicolas  Island  was  chosen.  This  study  uses 
data  frem  Vandenherg  AFB  since  it  is  a  coastal  station  and 
will  have  somewhat  different  surface  conditions  due  tc  the 
continental  influence.  See  Fig.  1.1  for  a  geographical 
orientation  of  the  area  of  interest.  Once  the  Vandenberg 
AFB  data  had  been  studied,  the  San  Nicolas  Island  results 
were  used  tc  make  comparisons.  The  results  of  this  compar¬ 
ison  are  of  interest  since  they  will  indicate  whether  or  not 
there  are  potential  problem  with  the  HABL  model  when  used  in 
a  regime  strongly  influenced  by  land  as  compared  to  a  regime 
closely  resembling  marine  conditions. 
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figure  1.1  Southern  California  Coastal  Region 


II.  HOD El  DESCBIPTIQN 


1.  E1CKGECQUD 

The  HAEL  model  is  use  at  the  Naval  Postgraduate  School 
(NPS)  is  a  zero-order,  two  layer,  integrated  mixed  layer 
model  (Davidson  et  al.  1984).  This  is  a  modified  version  of 
the  model  developed  by  Stage  and  Businger  (1981).  The  two 
layers  consist  of  the  wall-mixed,  turbulent  boundary  layer 
and  the  relatively  non- turbulent  free  atmosphere. 
Separating  these  layers  is  the  transition  zone  known  as  the 
inversion.  The  tars,  "zero-crder",  applies  to  this  model 
because  the  inversion  is  considered  a  jump  from  one  value  to 
another  instead  of  having  a  finite  gradient.  Thus,  we 
assume  that,  at  the  inversion,  the  equivalent  potential 
temperature  (9g),  and  the  total  specific  humidity  (q t  ) 
undergo  zerc-thickness  jumps. 

Total  specific  humidity  and  equivalent  potential  temper¬ 
ature  are  conserved  quantities  in  pseudo-adiabatic  processes 
and,  for  this  reason,  are  assumed  well-mixed  within  the 
boundary  layer.  They  ara  given  by: 

*t  =  qv  +  ql  (2.1) 

where  q  is  the  water  vapor  mixing  ratio  and  q  is  the 

*  L 

liquid  water  mixing  ratio;  and 

9=0  +  L/(Cq)  /I  21 

e  v  p  v  I 

where  9  is  the  potential  temperature,  Lv  is  the  latent  heat 
for  vaporization  of  water,  and  Cp  is  the  specific  heat  at 
constant  pressure  of  air.  A  third  fundamental  dependent 
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Figure  2. 1 


Typical  potential  teaperatnre  and 
moisture  profiles. 


variable  is  the  mean  mixed  layer  dapth  (h)  .  Typical 
profiles  of  these  parameters  are  shown  in  Fig.  2. 1.  As  can 
be  seer  in  Fig.  2.  1,  both  9e  and  qt  are  assumed  to  have 
constant  values  throughout  the  mixed  layer.  Above  this  is 
the  interfacial  regicn,  or  inversion,  where  the  jump  occurs. 

One  cf  the  more  appealing  features  about  this  model  is 
that  it  runs  from  input  parameters  gleaned  frcm  routine 
metecrclcgical  and  oceanographic  data.  These  parameters 
are:  surface  layer  humidity,  temperature  and  winds, 

vertical  gradients  cf  temperature  and  humidity  above  the 


mixed  layer,  ar.d  sea-surf acs  temperature  (SST)  .  F.eadily 
available  rawiascnde  observations  can  be  used  tc  detoiaiar 
ail  input  parameters  except  SST.  SST  values  car.  be  era  air.  id 
from  satellite  interpretation  cr  direcr  measurements. 

The  predictive  equations  for  the  three  dependent  vari¬ 
ables  are  given  by: 

dqt/dt  =  h  1(  (w'q;)Q  +  WgAqt) 

d6e/dt  =  h  (We;)0  +  Rc  -  RB  +  weA0e  )  (2.3) 

dh/dt  =  w  +  w 
D  e 

where  d/dt  represents  the  time  rate  of  change  of  the  vari¬ 
able  following  a  column  of  air.  (w'e')0  and  (w’q^Q  are  the 
surface  turbulent  fluxes  of  0  and  g  ,  while  the  A  terms  are 
the  magnitude  of  the  jumps  at  the  inversion.  Other  quanti¬ 
ties  introduced  here  are:  1)  we  -  the  entrainment  rate; 

2)  *D  -  the  large-scale  vertical  velocity  at  the  inversion 
height  (initialized  using  the  subsidence  rates  described 
later  in  this  chapter);  3)  P.B  -  the  rata  of  radiative  heat 
loss  per  unit  area  near  the  cloud  top;  and  4)  Ec  -  the 
rata  of  radiative  heat  gain  per  unit  area  at  the  cloud  base. 
(Haggerty,  1983) 

B.  NFS  MODEL  P&BTICQ1&BS 

Tie  use  of  a  well-mixed  layer,  along  with  other  develop¬ 
ments,  made  modifications  to  the  modal  described  by  Stage 
and  Eusinger  (1981)  necessary.  The  NPS  integrated  MABL 
differs  from  others  in  the  specifics  of  surface  flux,  radia¬ 
tive  flux  and  the  zero  gradient  assumption.  The  following 
discussions  are  brief  since  the  references  (Davidson  et  al. 
1984)  are  quire  complete. 
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1-  Jtl cs£heric  st-uctu^e 


The  profile  cf  a  well-mixed  atmospheric  variable,  $  , 
is  represented  as 


4>  ( 2 )  =  $ 

+  I\  2 

Z<  h 

ms 

*U)  =  tmh 

+  r. (2-h) 

9 

Z>  h 

(2.4) 


where  m  refers  to  the  mixed  layer,  s  to  the  surface  value,  h 
is  the  height  of  the  mixed  layer,  2  is  the  height,  A$  is  the 
jump  in  4>  across  the  inversion  zone,  is  the  gradient 
above  the  inversion,  and  r^m  is  the  gradient  in  the  mixed 
layer  (assumed  tc  be  zero)  . 

2.  Dynamics  within  the  mixed  layer 


The  model  predicts  the  time  evolution  of  the  mixed 
layer  values  of  g  ,  9  ,  ar.d  h  (see  Egn.  2.3).  At  initial!- 
zation,  the  lapse  rates  for  these  quantities  are  held 
constant  while  larger,  synoptic  scale  observations  provide 
the  input  values  of  the  boundary  layer  wind  for  the  entire 
forecast  period. 

Surface  fluxes  of  sensible  heat  and  water  vapor  are 
determined  through  an  application  of  the  bulk  aerodynamic 
formula  (Davidson  et  al.  1 S78) .  These  fluxes  are  driven  by 
the  bulk  air/sea  differences  and  are  expressed  as 


'"'Vo  ■  c5  l9o  -  V". 
'“'Vo  ■  c0  '*0  -  V". 


(2.5) 


where  Cg  is  the  ten  meter,  stability  dependent,  drag  coeffi¬ 
cient,  u*  is  the  friction  velccity,  and  9Q  and  qQ  are  the 
respective  sea-surface  values. 


rhar 


Ear amet erizaticn  for  the  entrainment  rate  is 
given  by  Stage  and  Eusinger  (1981)  .  This  paramete riza'rior. 
aliens  closure  of  the  system  cf  equations  and  determines  the 
time  evolution  of  the  inversion  height.  It  also  brings  the 

turbulent  kinetic  energy  (TKE)  budget  into  the  formulation 
since  entrainment  of  heat  flux  across  the  stable  inversion 
requires  TKE,  that  is,  work  is  done  by  this  process.  Here 
the  assumption  used  for  closure  states  that  the  dissipation 
rate  cf  TKE  is  a  fixed  fraction  (1-A)  of  the  production 
rate,  where  A  is  the  entrainmsnt  coefficient,  set  equal  to 
0.2.  Within  the  TKE  equation  are  terms  which  provide  for 
producticn  and  dissipation,  therefore  it  is  assumed  that 
both  cccur  withir  the  mixed  layer,  rather  than  cnly  where 
there  is  negative  buoyancy  flux.  These  are  associated 
through  the  use  of  the  turbulent  scaling  velocity  for  the 
mixed  layer,  w* ,  and  the  mixed  layer  depth  (h)  for  the 

scaling  length. 

3  .  gadiativa  flux  considerations 

when  considering  radiation  flux  and  its  effects  on 
the  mixed  layer  there  are  several  points  which  must  be 
understccd.  Distinction  must  be  made  between  leng-wave 

fluxes  for  clear  or  cloudy  cases,  and  short-wave  fluxes. 

a.  Long-wave  radiative  flux 

Onder  cloudless  sky  conditions,  the  net  long¬ 
wave  radiative  flux  is  determined  using  the  water  vapor  and 
temperature  profiles.  The  flux  is  calculated  by  inte¬ 
grating,  from  the  surface  (Z  =  0)  to  the  top  of  the  mixed 

layer  (Z*h)  ,  the  flux  emissivity  profile  where  the  emis- 

sivity  is  determined  from  the  Stef an-Boltzmann  law  (Fleagle 
and  Eusinger  1980) .  The  aforementioned  modification  has 
alsc  been  applied. 


The  presence  of  clcuds  necessitates  the  calcula¬ 
tion  of  the  net  radiative  flux  at  both  the  cloud  top  ar. a  the 
cloud  base.  These  calculations  are  also  made  using  the 
Stefan- Ecltzmann  law.  When  the  cloud  top  flux  is  being 
determined,  the  input  temperatures  are  the  cloud  top  temper¬ 
ature  and  the  effective  radiative  sky  temperature.  Per  the 
cloud  base  the  input  temperatures  are  the  temperatures  at 
the  cloud  base  and  at  the  sea  surface.  It  is  also  important 
to  note  that,  in  this  case,  the  flux  divergence  between  the 
sea  surface  and  the  cloud  base  is  neglected. 

t.  short-wave  radiative  flux 

The  H&BL  model  calculations  of  short-wave  radia¬ 
tion  take  into  consideration  both  direct  and  diffuse  compo¬ 
nents.  Eoth  scattering  and  absorption  contribute  to 
short-wave  extinction.  The  scattering  contribution  tc  the 
short-wave  radiative  flux  is  determined  through  use  cf  the 

delta-Eddington  method,  where  the  mixed  layer  scattering  is 
•  • 

due  tc  atmospheric  particles,  cloud  droplets  in  the  cloudy 
case  and  sea-salt  aerosols  in  the  clear  sky  case. 
Absorption  is  primarily  the  result  of  water  vapor  and  is 
taken  into  account  when  the  incident  flux  at  the  top  of  the 
mixed  layer  is  determined. 

C.  HCBEI  AEPLIC  AT 10 I 

All  model  runs  for  this  thesis  were  done  such  that  the 
situation  approximated  to  a  real  forecast  run.  That  is, 
the  mcdel  was  initially  set  tc  run  with  a  single  sounding 
and  an  estimate  of  -0.0005  m/s  for  the  subsidence.  Further 
runs  were  made  as  each  new  sounding  could  be  added.  This 
presents  a  realistic  representation  of  what  one  could  expect 
from  the  mcdel  for  the  Vandenberg  afb  area. 


Preparation  for  running  the  model  requires  that  several 
points  or  items  he  considered.  These  consist  of  ider.t i dying 
the  synoptic  pattern  and  initialization  parameters;  will  he 
discussed  below. 

1  *  Identif y  synoptic  Pattern 

Eefcre  checking  on  input  data  for  the  model  the 
overall  syncptic  pattern  must  be  checked.  To  do  this  cr.s 
utilizes  all  available  meteorological  data  (e.g,  surface  and 
upper  air  synoptic  analyses,  as  well  as  satellite,  rawin- 
sonde  and  individual  surface  observations).  The  desired 
synoptic  situation  fcr  optimum  model  predictability  is  to 
have  the  area  cf  interest  dominated  by  a  high  pressure 
syste i . 

2.  Identify  Initializa tion  Parameters 

The  next  step  in  preparing  to  make  the  model  runs  is 
to  deternine  the  input  parameters.  Potential  temperature 
and  moisture  values  fcr  the  mixed  layer,  height  and  strength 
cf  the  jump,  surface  wind,  sea-surface  temperature,  vertical 
velocity  and  the  lapse  rate  above  the  jump  are  the  necessary 
parameters  and,  of  these,  a  good  determination  of  the  jump 
value  across  the  inversion  is  critical. 

Except  for  vertical  velocities,  parameters  are 
determined  from  the  vertical  soundings,  and  can  be  done 
canually  cr  through  the  use  of  a  digitizing  program.  A 
sample  cf  a  digitized  ravinsonde  appears  in  Fig.  2.2. 
Extreme  care  must  be  used  sc  that  the  values  determined 
represent  the  true  scunding  or  serious  error  can  result. 

Th6  vertical  velocity  or,  in  our  case,  the  subsi¬ 
dence  rate  for  the  forecast  start  time  must  be  determined. 
There  are  several  means  for  calculating  vertical  velocity 
but  ncne,  tc  date,  provide  the  accuracy  desired.  with  this 
in  mind,  the  method  which  has  proved  to  be  most  accurate  was 


used.  This  method  is  the  Q-method  which  is  an  expression 
for  vertical  valccitj  derived  from  the  vertical  integration 
cf  the  moisture  budget  equation  and  is  given  by: 

A (q  +  l)  (dh/dt)  -  h  < d (q  +  l ) /dt )  +  (w'q1) 

w  s  —  11  ■  -  . . . . . . . .  (2.6) 

A  (q  +  l )  +  B(h/2) 

where  g  and  1  are  vapor  and  liquid  water  contents,  and  6, 
their  vertical  gradient  (Gleason,  1982).  The  computer 
program  used  to  calculate  the  vertical  velocity  uses  current 
data  as  well  as  data  within  the  previous  24  hours.  The  use 
cf  previous  rawinsonde  data  is  allowed  since  the  Q-method 
equation  includes  total  change  with  respect  to  time  (d  /dt) 
for  the  inversion  height  and  the  total  moisture.  The 
vertical  velocities  output  from  the  pregram  are  for  each 
data  period  input.  These  values  are  then  averaged  to  get 


figure  2.2 


Sample  of  digitized  rawinsonde 


the  rate  fcr  the  ertire  period  with  the  final  sutsicer.ce 
rate  calculated  teinc  the  one  used  in  the  model.  The  or.ly 
difficulty  encountered  with  this  method  was  when  the  subsi¬ 
dence  inversion  height  was  increasing  at  a  rate  sufficient 
to  yield  a  positive  average.  When  this  occurred  the  model 
input  rate  was  set  tc  zero.  Table  I  gives  the  subsidence 
rates  (m/s)  used  for  the  model  runs. 


TABLE  I 

Subsidence  Bates  for  Model  Buns 


CATE 

GMT 

SUESIDNECE 

DATE 

GMT 

SUBSIDENCE 

July 

25 

00 

-0.C05 

August 

29 

00 

-0.005 

12 

-0.C03948 

12 

-0 . 0  008  39  3 

26 

00 

-0.003359 

30 

00 

0.00 

12 

-0.0008265 

12 

-0.0002389 

27 

00 

-0.01489 

31 

00 

-0.0004126 

12 

-0. CO  1553 

12 

0.0 

28 

00 

-0.002856 

September 

1 

00 

0.0 

2 

00 

-0.005934 

12 

-0.001676 

vbei 

13 

00 

-  0 . 005 

3 

00 

-0.000082 

12 

-0.005813 

12 

0.0 

14 

00 

-O.C02856 

4 

00 

-0.0006649 

12 

0.0 

12 

-0.0002713 

5 

00 

-0.0005856 

12 

-0.001524 

pril 

12 

00 

-0.C05 

6 

00 

0.0 

12 

-0.002019 

14 

00 

O.C 

12 

-0.000365  1 

15 

00 

O.C 

12 

0.0 

16 

00 

-0.C02208 

12 

-0.002988 

23 


III.  S TNOPTIC  PATTERN  DISCUSSION 


Fcur  data  sets  were  used  for  evaluation  of  model 
validity.  These  sets  include  a  summer  period  (25  through  27 
July  1977),  a  winter  period  (13  and  14  December  1977),  a 
spring  period  (13  through  17  April  1981),  and  a  fall  period 
(29  Aug  through  5  September  1981).  Large  scale  analyses 
were  for  as  close  tc  1200  Greenwich  Mean  Tims  (GMT)  as 
possible.  In  most  cases  the  satellite  pictures  will  be  an 
1745  GMT  since  these  pictures  best  display  the  area  of 
interest . 

This  discussion  will  be  of  the  overall  weather  picture 
for  the  eastern  North  Pacific  Ocean  and  the  western  United 
States.  Included  will  be  conditions  over  the  area  as  a 
whole  and  unique  or  differing  conditions  at  Vandenberg  AF3 
itself.  Satellite  pictures,  surface  weather  analyses,  and 
upper  air  analyses  (5C0  millibars)  will  be  used  to  determine 
the  synoptic  situation  of  each  period. 

A.  SUMMER  CASE,  PERIOD  I  (25-27  JOLT  1977) 

The  surface  and  5C0  millibar  (mb)  synoptic  patterns  for 
this  period  appear  in  Figs.  3.2  through  3.4.  During  this 
period  the  coastal  region  was  dominated  by  a  high  pressure 
systex  ever  the  North  Pacific  Ocean  with  the  thermal 
troughing  extending  from  Mexico  northward  into  the  high 
desert  region  of  Ari2cr.a  and  southeastern  California.  Such 
a  syncptic  pattern  usually  produces  weak  onshore  flew  alcng 
the  southern  Califorria  coastal  regions.  The  upper  air 
analyses  (500  mb)  showed  no  strong  westerly  flew  over 
California  during  the  entire  period. 


Curing  this  period  the  Pacific  high  pressure  system 
moved  northeastward  and  on  27  July,  weakened  es  a 
stationary  front  developed  off  the  northwest  coast.  Also  of 
importance  on  27  July  was  the  development  cf  a  low  pressure 
center  west  of  Los  Angeles,  California  Even  with  this 
development  the  predominate  flow  remained  onshore  in  the 
Vanderberg  AFB  region. 

The  rawinsonde  data  revealed  a  subsidence  type  inversion 
throughout  the  period,  see  Fig.  3.1.  There  was  alec  a 
diurnal  variation  in  the  inversion  base.  Each  day  the  base 
ranged  from  the  surface  at  1200  GMT  to  a  height  above  650 
feet  at  C00C  GMT.  An  example  of  this  diurnal  change  also 
appears  in  Fig.  3.1. 

The  satellite  pictures  for  this  period  appear  in  Figs. 
3.5  through  3.7.  Orly  patchy  areas  of  lew  clouds  along  the 
California  coast  appear  on  the  satellite  pictures.  These 
areas  were  in  the  San  Francisco  bay  area  and  the  coastline 
south  cf  Lcs  Angeles.  The  only  ether  area  cf  low  clouds  was 
two  to  three  degrees  longitude  west  of  Vandenberg  AFB. 
There  was  no  significant  changes  to  the  coastal  locations  of 
the  lew  clouds  while  the  area  over  the  water  had  some  dissi¬ 
pation  and  was  even  further  offshore, 
vicinity  remained  under  clear  skies 


The  Vandenberg  AFB 
during  the  entire 


figure  3.5  Satellite  picture  for  25 


July 


30 


S.  11SI1B  CAS  E,  PEBICD  II  (13-14  DECEMBER  1977) 


The  general  surface  and  500  mb  analyses  appear  ir. 
figures  3.S  and  3.10  A  difference  between  this  synoptic 
pattern  and  that  found  in  the  other  periods  is  in  this 
patxern  the  desired  onshore  flow  exists  but  it  dees  not 
arise  from  the  strong  Pacific  high  and  thermal  troughirg  tut 
from  a  high  pressure  cell  that  had  separated  from  the  main 
Eacific  cell  and  moved  in  over  central  California  The 
strongest  westerly  flow  aloft  was  found  to  the  north  of 
Califcrria  with  the  Eclar  jet  moving  across  the  Washington  - 
Cregct  region.  Westerly  flow  over  Vandenberg  AFB  was 
stronger  than  in  other  periods. 

The  primary  aspect  is  that  on  13  December  a  high  pres¬ 
sure  center  moved  over  central  California  with  the  large 
high  pressure  system  still  maintained  over  the  North  Eacific 
Ccean.  Weak  thermal  troughing  was  also  extending  northward 
into  California  The  overall  result  of  these  surface 
patterns  was  weak  onshcre  flew  in  the  Vandenberg  AFB  and  San 
Nicolas  Island  regions.  The  only  other  significant  feature 
in  the  13  Eecemher  analysis  was  the  frontal  system  moving 
into  the  Pacific  Northwest. 

Cn  14  December  the  high  pressure  center  dissipated 
leaving  the  Eacific  high  pressure  system  and  ridging  acrcss 
central  California.  Weak  thermal  troughing  continued  in  the 
southern  California  area.  Offshore  flow  also  continued  to 
dominate  the  southern  California  coastline.  The  frcntal 
system  free  the  13th  moved  eastward  and  a  trailing  warm 
front  moved  to  northern  California. 

Fig.  3.8  is  an  example  of  the  vertical  soundings  for 
this  period.  The  rawinsonde  data  showed  subsidence  type 
inversions  on  both  days.  Most  important  is  that  on  13 
December  at  0000  GMT  the  inversion  was  very  weak  and  nox 
well  structured.  The  inversion  strengthened  and  remaired 
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aloft  until  1200  GMT  cn  14  December  when  it  was  observe!  as 
a  surface  based  inversion.  At  this  time  the  inversion 
reached  the  maximum  strength  for  the  period. 

The  satellite  pictures  for  this  period  appear  ir.  Figs. 
3.11  and  3.12.  Clear  skies  in  the  Vandenberg  A FB  area  with 
a  frontal  system  moving  intc  the  northern  California  regions 
appear  in  the  satellite  picture  for  13  December.  There  is 
very  little  indication  of  low  clouds  along  the  southern 
California  coastline.  By  14  December  the  front  is  well  into 
northern  California  with  increased  cloudiness  across  central 
California.  There  has  also  been  an  increase  in  the  amount 
of  lew  clouds  along  the  southern  California  coast.  patchy 
stratus  fellows  the  coast  south  cf  Los  Angeles  and  a  large 
area  now  sits  off  the  Vandenberg  AFB  coast. 


figure  3.8  Example  cf  Sinter  Case  rawinsonde 


figure  3.10  Surface  (a 
1200  GIT  14 


Figure  3.11 


Satellite  picture  for  13  December  1977. 
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1745  14BE77  32 A- 1  00231  16051  WB2 


Figure  3.12  Satellite  picture  for  14  December  1977 
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C.  SEEING  CAS E#  PEBICD  III  (13-16  APHIL  1981) 


The  synoptic  patterns  for  this  period  appear  in  Figs. 
3.14  through  3.17.  The  dominating  feature  was  a  high  pres¬ 
sure  system  extending  northeastward  from  the  North  Pacific 
Ocean  regions  across  Washington  into  Canada.  Thermal 
troughing  dominated  the  high  desert  regions  of  California 
and  extended  into  northern  California.  The  predominate 
southern  California  coastal  flow,  generated  by  these 
patterns,  was  onshore.  500  mb  analyses  shew  light  and  vari¬ 
able  winds  with  tendencies  toward  formation  of  a  lew  pres¬ 
sure  center  through  the  period.  The  strongest  westerly  flow 
remained  well  north  of  the  Var.denberg  AFB  area. 

Cn  13  April  the  Eacific  high  was  strong  but  showing  a 
tendency  to  weaken  as  the  longwave  pattern  extended  into 
Canada.  As  the  longwave  high  pattern  shifted  to  the  central 
United  States  the  Pacific  high  weakened.  The  thermal  trough 
in  southern  California  continued  to  dominate  but  also  had 
some  weakening. 

Ey  15  April  the  Pacific  high  system  had  begun  to 
strengthen  with  some  ridging  across  central  California.  The 
high  continued  intensification  through  the  remainder  of  the 
period.  Cn  16  April  the  ridging  extended  into  the 
Washington  area  with  the  Pacific  high  pressure  center  begin¬ 
ning  tc  move  northeastward.  Throughout  this  portion  of  the 
period  the  thermal  trough  showed  very  little  variation.  The 
only  ether  key  occurrence  from  the  15th  to  the  16th  was  the 
movement  of  a  cold  front  across  the  Pacific  Northwest  into 
Nevada. 

Fig.  3.13  is  an  example  of  the  rawinsonde  soundings  for 
this  period.  Hawinsonde  reports  during  this  period  had 
subsidence  inversions  with  much  daily  variation.  At  the 
beginning  of  the  period  the  inversion  was  weak  and  aloft. 
Eiurnal  variations  intensified  the  inversion  but  also 
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brought  it  tc  the  surface.  This  pattern  continued  through 
the  14th  but  on  the  15th  the  inversion  intensified  at  J COO 
GMT  and  remained  alcft  until  1200  GMT  on  the  16th.  Durir.o 
the  latter  half  cf  the  period  there  was  some  fluctuation  in 
the  intensity  but  over  all  the  inversion  remained  strong. 

These  satellite  pictures  appear  in  Figs.  3.18  through 
3.21.  On  13  April  the  satellite  showed  Vanderberg  APE  under 
clear  skies  while  the  coastal  region  south  of  there  was 
dominated  by  low  clouds.  Ey  14  April  the  area  of  coverage 
fcy  low  clouds  increased  and  Vandenberg  AFB  was  within  that 
area.  low  clouds  continued  to  dominate  the  Vandenfcerg  AFB 
area  with  seme  dissipation  beginning  on  15  April.  16  Apr 
showed  crly  thin  coverage  of  low  clouds  over  the  Vandenberg 
AFB  area  with  the  denser  coverage  south  of  the  area. 


SPEC.  MUM.  Igskg)  5 
POT.  TEMP.  (C)  0 


Figure  3.13  Example  of  rawinsonde  for  Spring  Case. 


figure  3.14 


figure  3.19  Satellite  picture  for  14  April  1981. 
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1745  15AP31  36A-2  00243  16002  WB2 


figure  3.20  Satellite  picture  for  15  April  1981. 
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1745  1 SAPS 1  36A-2 


00251  15382  WB 
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0.  I  ILL  CASE,  P  SB  IOC  IV  (29  AUGUST  -  5  SEPTEMBER  1981) 

Ihe  surface  (a)  and  500  at  (b)  analyses  for  Period  IV 
appear  in  Figs.  3.23  through  3.30.  The  predominate  surface 
feature  during  this  extire  period  was  rhe  presence  of  a  high 
pressure  system  in  the  Eastern  portions  of  the  North  Pacific 
Cesar.  A  thermal  low  and  trough  system  dominated  the 
southern  California  region.  The  effect  of  these  two  surface 
features  was  to  produce  weak  onshore  flow  in  the  area  of 
concern.  Conclusions  drawn  from  the  500  mb  analyses  show 
that,  except  for  the  last  day  of  the  period,  the  maximum 
westerly  flow  remained  north  of  California.  Only  cn  5 
September  did  the  Pclar  jet  extend  toward  the  regions  of 
interest . 

The  surface  analysis  for  29  August  (Fig  3.23)  shews  the 
western  Pacific  Ocean  dominated  by  a  high  pressure  system 
with  a  lew  pressure  center  located  off  the  Washington  ccast. 
A  thermal  low  extended  northwestward  from  Mexico  into 
Calif erria  with  troughing  extending  into  northern 
California.  The  southern  end  of  a  cold  front  had  also  moved 
into  the  northern  California  area. 

Ey  30  August  the  cold  front  had  moved  into  southern 
Califcrria  and  the  thermal  low  had  changed  its  orientation. 
The  Facific  high  pressure  system  showed  weak  intensification 
with  movement  toward  the  northeast  and  had  pushed  inland  by 
the  31st.  The  31st  also  shewed  that  the  cold  front  was 
beginning  to  leave  California  but  with  another  frontal 
system  moving  out  of  the  Gulf  of  Alaska.  Figures  3.24 
through  2.25  depict  these  occurrences. 

The  1  September  analysis  revealed  a  weakening  of  the 
high  pressure  and  a  cold  front  approaching  California.  A 
thermal  trough  was  well  established  over  California  tut  very 
diffuse  cn  2  September  as  the  cold  front  moved  through 


southern  California.  The  high  pressure  system  again  inten¬ 
sified  ar.d  showed  r.o  marked  weakening  throughout 
remainder  cf  the  period.  The  only  significant  event  after  2 
September  was  the  entrance  of  a  rapidly  moving  cold  front 
into  California.  This  scenario  is  shown  in  Figs  3.26 
through  3. 30. 

A  saiple  of  the  rawinsonaes  for  this  period  appears  in 
Fig.  3.22.  The  rawinscndss  for  this  period  have  a  very  good 
subsidence  type  inversion  throughout  the  period.  In  each 
rawinscrde  release  the  inversion  was  found  above  1000  feet 
with  decreasing  relative  humidity  through  the  inversion 
zone. 

The  extent  and  development  of  the  low  clouds  appear  in 
the  pictures  in  Figs.  3.31  and  3.38.  The  satellite  pictures 
revealed  a  buildup  cf  coastal  low  clouds  (stratus  and/or 
fog)  during  the  period.  Beginning  on  29  August  the  clouds 
tended  to  be  thin  with  seme  dissipation  during  the  day. 
Further  intc  the  period  the  stratus  thickens  and  by  the  end 
of  the  period  there  is  a  dense  area  of  low  clouds  off  the 
central  and  southern  California  coastal  regions. 
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figure  3.23  Serf  ace  (a)  and  500  ab  (b)  analyses 

1200  GBT  29  lugust  1977. 


figure  3.24  Surface 
1200  GHT 


% 


) 

and  500  eb  (b)  analyses 
lugust  1977. 


(b> 

figure  3.26  Surface  (a)  and  500  mb  (b)  analyse 
1200  GBT  1  September  1977. 


(b) 

Figure  3.29  Surface  (a)  and  500  ab  (b)  analyses 
1200  GHT  *  September  1977. 

58 


M* 


002 x 3  15992  U82 


figure  3.33  Satellite  picture  for  31  August. 
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Figure  3.34  satellite  picture  for  1  September  1981. 


2  16011  WB2 


5  Q3SE31  36A-2 _ DO 


Figure  3.36  Satellite  picture  for  3  September  1981. 
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Figure  3.38  Satellite  picture  for  5  September  1981 
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IV.  BESOLTS 


Bedel  forecasts  were  made  for  every  30  minutes  for  24 
hours  following  the  initialization  time.  The  forecasts 
presented  in  the  following  tables  show  only  the  12  and  24 
hour  values  since  they  correspond  no  the  verification  times. 
The  structure  of  the  tables  are  such  that  the  first  fore¬ 
cast,  indicated  by  column  beading  FCST1,  was  based  on  the 
first  observation  for  initialization  and  within  that  column 
are  the  forecasts  generated.  For  example,  if  the  first 
observation  was  for  1/0000  GMT  then  the  forecasts  from  that 
observation  would  be  for  1/1200  GMT  and  2/0000  GMT.  These 
forecasts  would  be  in  column  FCST1.  The  second  forecast, 
column  beading  FCST2,  was  initialized  with  the  second  obser¬ 
vation  (e.g.,  1/1200  GMT)  and  within  that  column  are  the 
forecasts  generated  (2/0000  GMT  and  2/1200  GMT).  Each  table 
continues  it  this  manner  until  all  forecasts  for  a  period 
have  beer  given. 

Boot  mean  square  (RMS)  scores  are  given  following  each- 
forecast  except  the  last.  No  RMS  for  the  last  forecast  were 
given  as  there  were  no  observations  for  comparison.  Also 
note  that  the  12  heur  RMS  value  simply  represents  the 
difference  between  observed  and  forecasted  values. 

The  values  given  in  the  tables  are  for  the  parameters: 
inversion  height,  lifting  condensation  level,  mixed  layer 
potential  temperature  and  mixed  layer  specific  humidity. 
Dnits  for  these  are  meters  for  the  inversion  height  and  the 
lifting  condensation  level,  degrees  Centigrade  for  potential 
temperature,  and  g/kc  for  specific  humidity.  Also  nets  that 
subsidences,  used  in  discussing  these  results,  are  given  in 
Table  I. 
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a.  DISCUSSION  OF  MODEL  RESGLTS  FOR  PERIOD  I  (SOMBER) 

This  period  was  dcminated  by  clear  skies.  The  forecast 
lifting  ccndensa  tion  level  values  (see  Table  III)  cr.ly  fore¬ 
casted  clcuds  during  PCS  12.  The  forecasted  inversion 
heights.  Table  II,  reveal  that  the  model  was  either  very 
close  or  quite  in  error  relative  to  the  observed  height. 
Tables  I1?  and  V  reveal  that,  throughout  the  period,  there 
were  alsc  problems  associated  with  potential  temperature  and 


T ABLE  II 

July  Inversion  Results 


DATE 

GET 

CBS 

POST  1 

RMS 

FCST2 

RMS 

FCST3 

RMS 

25 

00 

356 

— 

- 

— 

— 

— 

- 

12 

157 

195 

38.0 

- 

- 

- 

- 

26 

00 

277 

108 

173.2 

92 

185.0 

- 

- 

12 

199 

— 

- 

112 

204.4 

165 

34.0 

27 

00 

206 

— 

— 

— 

— 

100 

111.3 

DATE 

GMT 

CBS 

FCST  4 

RMS 

FCST5 

RMS 

FCST6 

RMS 

FCST7 

21 

CO 

206 

172 

34.0 

— 

— 

— 

12 

120 

146 

42.8 

154 

34.0 

— 

- 

- 

28 

cc 

199 

— 

— 

116 

89.7 

71 

128.0 

- 

12 

— 

— 

- 

- 

— 

42 

1  10 

29 

CO 

- 

- 

■  - 

- 

- 

- 

- 

62 

TAELE  III 

July  Lifting  Condensation  Levels 


DATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

25 

00 

633 

- 

- 

— 

• 

12 

146 

2  99 

— 

— 

- 

26 

00 

1119 

199 

42 

— 

- 

12 

217 

- 

0 

519 

- 

27 

00 

206 

— 

— 

309 

185 

12 

452 

- 

- 

— 

168 

DATE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

57 

12 

452 

605 

— 

— 

28 

00 

1261 

410 

184 

- 

12 

- 

- 

176 

586 

29 

00 

- 

- 

- 

366 
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T  ABLE  I? 


July  Potential  Temperatures 


DATE 

GMT 

OBS 

FCST1 

FC5T2 

FCST3 

FCST4 

25 

00 

16.2 

- 

— 

— 

• 

12 

11.  6 

17.  a 

— 

- 

— 

26 

00 

17.4 

13.2 

11.6 

- 

- 

12 

12.  6 

• 

8.5 

17.8 

— 

21 

00 

20.9 

- 

— 

18.4 

15.  2 

12 

14.  1 

— 

- 

- 

16.7 

DATE 

GiiT 

OBS 

FCST5 

FCSTo 

FCSI7 

27 

12 

14.  1 

19.8 

— 

- 

23 

00 

19.5 

19.6 

17.1 

- 

12 

— 

- 

18.3 

18.9 

29 

00 

- 

- 

- 

19.3 

TABLE  ? 

July  Boisture  Values 


DATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

25 

00 

8.2 

- 

- 

* 

12 

7.9 

6.3 

— 

— 

- 

26 

00 

6.8 

5.4 

8.3 

- 

— 

12 

8.  1 

- 

7.2 

9.7 

• 

27 

00 

9.3 

- 

11.3 

9.8 

12 

7.  9 

- 

— 

- 

10.8 

BATE 

GMT 

OBS 

FCST5 

ECST6 

FCST7 

27 

12 

7.9 

10.6 

1 

28 

00 

7.3 

11.6 

11.  1 

12 

— 

- 

12.0 

10.1 

29 

00 

- 

- 

- 

11.6 

specific  humidity  forecasts.  These  problems  will  be  covered 
during  the  following  discussions. 

Examination  of  the  inversion  height  results  reveals  seme 
good  and  seme  poor  height  forecasts  at  the  12  hour  point  but 
mostly  bad  forecasts  at  24  hours.  Those  forecasts  made  with 
average  subsidence  rates  (-  0.003  to  -0.005  m/s)  had  fair  to 
good  results  at  12  hours  but  had  poor  results  at  24  hours. 
The  best  height  forecasts  (Table  II  column  FCST4)  were  with 
the  weakest  subsidence  rate  while  the  second  best  (the  fifth 
forecast)  was  with  the  strongest  subsidence  rate.  Also  of 
interest  is  that  forecasts  with  the  strongest  subsidence 


10 


above 


also  had  the  liftirg  condensation  level  well  above  -he 
inversion.  Another  important  feature  in  Table  II  is  th:t, 
except  fcr  the  forecast  with  the  weakest  subsidence,  the 
test  forecasts  were  for  the  12  hour  periods  initialized  at 
1200  CM  1  (a  morning  period  fcr  the  California  region)  . 

Ectsrtial  temperature,  specific  humidity,  radiation 
(with  and  without  clouds)  and  entrainment  will  be  examined 
in  explaining  what  was  occurring  with  the  predicted  inver¬ 
sion  heights  for  this  period.  Each  parameter  will  be 
described  separately  as  to  its  possible  impact  on  the  inver¬ 
sion  forecast. 


1)  potential  temperature  -  Under  clear  skies  one 
would  normally  expect  nightly  cooling  ever  land. 
The  observed  values  shown  in  Table  IV  reveal  that 
this  occurred  yet  the  forecasts  do  net  reflect  this 
fact.  Nocturnal  ceding  strengthens  the  jump  and 
when  coupled  with  subsidence  lowers  the  inversion  as 
seen  in  Fig.  4.1.  All  forecasts  were  lowering  the 
inversion  due  to  good  subsidence  rates  for  initiali¬ 
zation.  Daytime  heating  acts  against  the  subsi¬ 
dence,  weakens  the  jump,  and  results  in  lifting  of 
the  inversion.  This  diurnal  event  was  not  reflected 
in  the  forecasts. 

2)  specific  humidity  -  After  the  first  two  fore¬ 
casts,  specific  humidity  forecast  values  were  only 
fair.  There  were  nc  forecast  values  which  were  good 
since  the  errers  ranged  from  just  under  1  g/kg  up  to 
3.3  g/kg.  The  primary  trend  in  the  forecasts  was  to 
increase  specific  humidity.  This  prediction  cf  too 
large  a  specific  humidity  caused  the  model  to 
predict  even  lower  lifting  condensation  levels 
which,  in  turn,  enhances  the  poor  inversion 
predictions. 
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3)  radiation  -  This  factor  aus-  be  considered  under 
conditions  of  clear  as  wall  as  cloudy  skies. 

a)  cloudy  -  In  the  second  forecast  the  sole! 
predicted  clouds  to  form  and  extend  tc  the 
surface.  The  clouds  radiate  long  wave  energy  to 
space  so  the  potential  temperatures  were  forecast 
too  low.  Inversion  predictions,  while  decreasing 
then  increasing  during  this  forecast,  did  not 
reflect  the  diurnal  change  predominate  with  clear 
skies . 

t)  clear  -  Radiation  computations  here  could 
have  been  incorrectly  influenced  by  the  large 
specific  humidity  forecasts.  This  accompanied 
with  good  subsidence  rates  caused  all  forecasts 
to  lower  the  inversion. 

4)  entrainment  rates  -  these  influence  the  time 
rate  of  change  of  specific  humidity  and  potential 
temperature  as  well  as  inversion  height.  In  several 
cf  the  forecasts  the  entrainment  rate  was  zero  so  it 
was  internally  set  to  a  positive,  yet  minimal, 
value.  This  could  have  affected  the  other  fore¬ 
casted  values  but  it  is  believed  the  effect  was  less 
than  others. 

ft  definite  diurnal  pattern  can  be  seen  (see  Fig.  4.1)  in 
the  observed  inversion  heights.  Generally  there  were 
increased  heights  in  the  transition  from  morning  to  evening 
and  decreases  fcr  the  opposite  transition.  This  pattern 
coincides  with  the  diurnal  pattern  of  the  potential  tempera¬ 
tures  which  reflect  radiative  daytime  heating  and  nightly 
cooling.  Surface  winds  also  displayed  the  influence  of 
daily  heating  with  evening  reports  having  the  strongest  wind 
speeds  and  with  steady  northwesterly  directions  (an  onshore 
wind).  Figure  4.1  also  depicts  the  12  and  24  hour  forecasts 
and  it  is  easily  seen  than  the  diurnal  pattern  shown  by  the 
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observed  data  is  net  shown  fcy  the  forecasts.  Only  tn-= 
nornirg  tc  evening  12  hour  forecast  indicate  a  correct 
trend.  However  close  inspection  of  the  forecasts  indicate 
that  at  thesa  times  the  model  inputs  for  subsidence  were 
strong  thereby  influencing  the  forecast. 

These  effects  all  suggest  problems  with  the  radiation 
calculations  throughout  this  summer  period.  A  reason  for 
this  error  is  the  fact  that  Vandenberg  AFB  is  a  shoreline 
station  and  the  model  is  primarily  designed  for  over  water 
use.  Therefore,  land  heating  and  cooling  by  solar  radia¬ 
tion,  even  through  a  sxratus  sky,  require  adjustments  to 
apply  the  model  at  a  shoreline  station.  Vertical  velocity 
and  entrainment  rates  also  affected  forecasted  MABL  proper¬ 
ties  through  influences  on  the  rate  of  change  of  key  parame¬ 
ters  as  well  as  the  inversion  heights  (see  Egn.  2.3). 

B.  DISCUSSION  OF  HOTEL  RESULTS  FOR  PERIOD  II  (WINTER) 

Results  in  Table  VI  generally  suggest  that'  during  this 
season  the  model  fails  to  forecast  variations  in  the  height 
cf  the  inversion.  Tables  VII  through  IX  compare  forecast 
and  observed  lifting  condensation  levels,  potential  tempera¬ 
tures  and  specific  humidity  values. 


TABLE  ¥1 

December  Inversion  Results 


DATE 

GBT 

CBS 

FCST  1 

RMS 

FCST2 

SMS 

FCST3 

RMS 

FCST4 

1  3 

CO 

164 

— 

— 

- 

- 

- 

— 

• 

12 

171 

45 

126.0 

- 

— 

- 

- 

— 

14 

00 

284 

14 

210.7 

40 

244.0 

- 

— 

- 

12 

297 

— 

— 

11 

265.8 

187 

110.0 

— 

15 

00 

— 

— 

— 

- 

- 

124 

492 

12 

- 

- 

- 

- 

- 

- 

- 

623 

7  4 


TABLE  VII 

December  lifting  Condensation  Levels 


DATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCS  TU 

13 

00 

1083 

- 

— 

— 

— 

12 

2055 

123 

— 

- 

— 

14 

00 

1056 

22 

936 

— 

- 

12 

174 

- 

338 

511 

— 

15 

00 

- 

- 

- 

255 

0 

12 

- 

- 

- 

- 

0 

TABLE  VIII 

December  Potential  Temperatures 


CATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

13 

00 

17.2 

— 

— 

— 

— 

12 

15.  8 

14*6 

- 

- 

— 

14 

00 

15.7 

14.  1 

15.4 

- 

- 

12 

13.0 

- 

15.2 

14.9 

— 

15 

00 

— 

— 

— 

14.8 

8.2 

12 

- 

- 

- 

- 

5.4 

I 


TABLE  IX 

December  Moisture 

Values 

CATE 

GMT 

OBS 

FCST1  FCST2 

FCST3 

FCST4 

13 

00 

8.9 

-  - 

- 

- 

12 

3.6 

9.7 

• 

- 

14 

00 

6.3 

9.9  6.6 

— 

— 

12 

8.  3 

9.0 

8.0 

— 

15 

00 

— 

-  — 

9.2 

7.3 

12 

- 

•  ^ 

- 

7.0 

Vandenfcerg  AFB  remained  under  clear  skies  in  this 
yet  the  forecasts  sere  for  rapid  formation  of 
Forecast  specific  humidity  values  improved  in  the 
forecast  hut  were  still  not  acceptable  for  the 
period.  The  potential  temperature  forecasts  were 
the  12  hour  point. 


period 
clouds, 
second 
ent  ire 
best  at 


The  primary  causes  for  the  poor  performance  of  the  modal 
in  this  case  appear  tc  be  vertical  velocity  and  entrainment, 
lhe  radiation  and  cloud  calculations  contributed  tc  thr- 
error,  tut  to  a  lesser  extent.  Entrainment  was  a  key  factor 
during  the  first  three  forecasts  because  null  values  were 
obtained  in  computation,  so  siniaal  positive  values  were 
assigted.  The  radiation  segment  calculations  are  strcnaly 
influenced  by  the  presence  cf  clouds  and  on  the  last  day  the 
modal  forecast  clouds.  With  zero  vertical  velocity,  the 
model  forecast  rapid  lifting  of  the  inversion.  With  no 
observations  it  is  hard  to  evaluate  those  forecasts  but  the 
trend  was  for  increasing  heights  although  not  as  rapidly  as 
predicted. 

forecast  and  observed  heights  for  this  period  appear  in 
Figure  4.2.  The  diurnal  variations  noted  in  the  summer  case 
are  net  present  in  this  case.  It  is  believed  that  the  large 
errors  in  the  inversion  height  forecasts  are  due  to  two 
factors,  first,  the  vertical  velocities  with  which  the  runs 
were  initialized  were  large  which  lead  to  forecasts  of 
lowering  heights  in  all  but  the  final  run.  In  reality,  the 
heights  shewed  a  continual  increase.  The  vertical  velocity 
used  in  the  last  run  was  such  that  rapid  lifting  of  the 
inversion  was  forecast.  Second,  the  surface  winds  which 
stayed  strong  even  through  the  morning  hours  allowed  mere 
mixing  in  the  lower  layers  thereby  increasing  the  depth  of 
the  fixed  layer.  This  is  especially  noticable  in  the 
13/1200  GMT  surface  wind  (360  degrees  at  2.1  m/s). 

At  Vandenberg  AFE  a  small  chain  of  mountains  exists 
due-north  of  the  base.  With  this  northerly  flow,  there 
would  be  an  increase  in  the  local  turbulence  in  the  lower 
levels.  Therefore,  the  poor  height  forecasts  in  this  case 
is  associated  with  the  land  influence  which  is  net  taken 
into  account  by  the  ircdel. 
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C.  DISCUSSION  OP  MODEL  HESOLTS  POB  PERIOD  III  (SPRING) 


Tables  X  through  XIII  list  the  forecasxed  values  for  the 
spring  case.  This  period  was  cloud  free  the  first  d 

had  lew  clouds  the  second.  The  forecasts  die  net  include 
cloud  formation  until  the  third  day.  Even  though  the  clcud 
forecasts  were  in  error  it  dees  not  appear  that  the  cloud 
calculations  had  any  real  effect  on  the  poor  inversion 
height  predictions. 

The  vertical  velocities  and  entrainment  rates  were 


TABLE  X 

April  Inversion  Results 


CATE 

GMT 

OBS 

FCST1 

RMS 

FCST2 

RMS 

FCST3 

RMS 

13 

00 

141 

— 

- 

— 

— 

— 

— 

12 

136 

31 

105.0 

- 

- 

— 

— 

14 

00 

289 

9 

299.0 

75 

214.0 

- 

— 

12 

150 

- 

40 

240.6 

293 

143.0 

15 

00 

469 

- 

- 

- 

- 

298 

222.9 

CATE 

GMT 

OBS 

FCST4 

RMS 

FCST5 

RMS 

FCST6 

RMS 

15 

00 

469 

140 

329.0 

— 

12 

509 

128 

503.4 

624 

115.0 

- 

— 

16 

00 

421 

835 

429.7 

6  84 

263.0 

/ 

12 

291 

- 

— 

- 

— 

867 

663.2 

CATE 

GMT 

OBS 

FCST7 

‘  RMS 

FCST8 

16 

12 

291 

340 

4  9,0 

17 

00 

— 

440 

- 

323 

12 

— 

- 

- 

335 

different  than  those  for  other  period.  Entrainment  rates 
were  zero  fer  the  first  and.  third  forecasts,  and  the  first 
three  hours  of  the  fifth  forecast  and  were  assigned  a 
nininal  positive  value  of  .001.  Hence,  the  inversion  height 
was  sutsidence  controlled  for  these  times.  The  moisture 
forecasts  were  good  since  values  agree  to  within  2  g/kg  of 
the  observed  value  and  in  some  cases  even  less  than  1  g/kg. 
Each  of  these  forecasts  were  initialized  with  the  evening 
soundirg  in  which  lifting  condensation  levels  indicated  no 
clouds. 


TiBLE  XI 

April  Lifting  Condensation  Levels 


EATE 

GXT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

13 

00 

705 

- 

— 

- 

- 

12 

3  70 

103 

- 

— 

— 

14 

00 

989 

59 

274 

— 

— 

12 

305 

— 

191 

919 

— 

15 

00 

506 

- 

- 

861 

234 

12 

445 

— 

• 

— 

190 

E  ATE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

FCST8 

15 

12 

445 

315 

• 

— 

- 

16 

00 

623 

423 

398 

- 

- 

12 

256 

— 

303 

297 

- 

17 

00 

- 

— 

— 

206 

78 

12 

- 

- 

- 

- 

0 

TABLE 

XII 

April 

Potential  Temperatures 

DATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

13 

00 

13.6 

— 

— 

— 

— 

12 

7.3 

12.8 

— 

— 

— 

14 

00 

15.7 

12.9 

10.7 

- 

- 

12 

8.2 

— 

12.5 

15.1 

- 

15 

00 

13.  2 

• 

« m 

- 

14.7 

10.3 

12 

10.6 

— 

— 

11.4 

EATE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

PCS  T8 

15 

12 

10.6 

9.9 

- 

— 

- 

16 

00 

12.  4 

9.6 

8.7 

- 

— 

12 

8.3 

— 

7.  1 

12.3 

— 

17 

00 

— 

— 

— 

9.8 

5.6 

12 

- 

- 

- 

- 

2.7 

Close  examination  reveals  that  the  potential  temperature 
forecasts  did  not  reflect  the  observed  diurnal  variaticn. 
The  12-hcur  forecast  generally  showed  cooling  or  heatir.q  as 
it  shculd  but  not  sufficiently  to  agree  with  observed 
results.  The  24-hour  forecast  tended  to  reverse  the  diurnal 
trend  and  to  have  coding  instead  of  heating  cr  vice-versa. 

It  is  believed  that  the  land  again  influenced  the 
observed  inversion  heights.  Figure  4.3  shows  the  observed 
and  forecasted  inversion  height  values  for  this  case.  Close 
examination  reveals  that  there  appears  to  be  two  problems 
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TABLE  XIII 

April  Hoist ure  7alaes 


CATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

13 

00 

6.6 

— 

— 

- 

• 

12 

5.2 

8.9 

— 

- 

- 

14 

00 

6.5 

9.0 

6.9 

- 

— 

12 

5.7 

- 

8.2 

6.5 

- 

15 

00 

7.2 

— 

- 

6.6 

6.8 

12 

6.2 

— 

- 

— 

7,6 

EATE 

GMT 

OBS 

FCST5 

ECST6 

FCST7 

FCST8 

15 

12 

6.2 

6.4 

— 

- 

— 

16 

00 

6.4 

5.9 

5.6 

— 

- 

12 

5.9 

— 

5.3 

7.6 

- 

17 

00 

— 

— 

- 

6.7 

5.4 

12 

- 

- 

- 

- 

5.2 

which  are  a  combination  of  the  problems  discussed  for  the 
summer  and  winter  cases.  The  first  part  of  the  period  was 
strongly  influenced  fcy  diurnal  fluctuations  and  the  fore¬ 
casted  values  were  similar  tc  those  in  the  summer  case.  The 
winter  case  was  influenced  by  subsidence  and  surface  winds. 
In  this  case,  it  isr  believed  the  northerly  winds  caused  the 
increase  in  the  observed  inversion  at  15/0000  GMT  due  to 
tcpcgraphic  induced  large  mechanical  mixing.  The  vertical 
velocity  allowed  the  model  to  decrease  the  heights  for  that 
time . 

I.  El 'CESSION  OF  HOE El  BESOLTS  FOB  PEBIOD  IV  (FALL) 

Tables  XIV  -  XVII  list  the  observed  and  forecasted 
values.  As  before  predicted  inversion  heights  for  this  case 
appear  tc  be  no  better  than  those  for  previous  cases.  The 
formation  of  clouds  was  forecasted  accurately  since  the 
modal  predicted  cloud  formation  at  a  time  when  low  clouds 
cccurred  in  the  area.  However,  forecasted  lifting  condensa¬ 
tion  levels  were  not  in  agreement  with  those  observed. 
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Aug ust-September  Inversion  Results 


CATE 

GMT 

OBS 

FCST1 

RMS 

FCST2 

RMS 

FCST3 

RMS 

29 

CO 

287 

- 

- 

— 

- 

- 

- 

12 

297 

137 

160.0 

— 

- 

- 

- 

30 

00 

365 

67 

336.2 

524 

159.0 

- 

— 

12 

31  1 

- 

— 

6  38 

363.6 

3  83 

72.0 

21 

CO 

391 

— 

- 

— 

- 

4  96 

127.3 

CATE 

GMT 

OBS 

FCST4 

RMS 

FCST5 

RMS 

FCST6 

RMS 

2  1 

00 

39  1 

512 

12  1.0 

— 

12 

553 

666 

165.6 

574 

21.0 

- 

- 

1 

00 

322 

- 

— 

755 

433.5 

819 

497.0 

12 

M 

— 

— 

— 

— 

1057 

- 

CATE 

GMT 

OBS 

FCST7 

RMS 

FCST8 

RMS 

FCST9 

RMS 

1 

12 

a 

363 

— 

— 

— 

— 

2 

00 

356 

£45 

18S.0 

— 

— 

— 

- 

12 

356 

- 

- 

273 

83.0 

- 

- 

3 

00 

4  1  1 

— 

- 

225 

244.2 

4  49 

38.0 

12 

469 

— 

— 

- 

5  09 

55.2 

CATE 

GMT 

OBS 

FCST10 

1  RMS 

FCST11 

RMS 

FCST12 

RMS 

3 

12 

469 

491 

22.0 

— 

— 

- 

4 

00 

422 

714 

292.8 

684 

262.0 

— 

— 

12 

529 

- 

- 

8  54 

417.5 

454 

75.0 

5 

00 

424 

— 

— 

— 

— 

562 

157.1 

CATE 

GHT 

OBS 

FCST13 

RMS 

FCST14 

RMS 

FCST15 

5 

00 

424 

675 

2510 

— 

12 

429 

823 

467. 

2  404 

25 

- 

6 

00 

— 

- 

— 

496 

4  98 

12 

- 

- 

- 

561 

Ihe  forecasted  specific  humidity  values  have  an  inter¬ 
esting  trend  (shewn  in  Table  XVII)  .  Forecasts  initialized 
from  evening  data  were  good  but  those  initialized  from 
morning  data  were  bad.  When  initialized  with  evening  data, 
the  predicted  specific  humidity  values  were  representative 
cf  the  ctserved  values.  Obviously,  pcor  forecasts  of 
specific  hunidity  or  any  other  key  parameter  lead  to  pcor 
prediction  of  inversion  heights. 

Ectertial  temperatures,  in  all  but  two  runs,  were  fore¬ 
cast  to  always  decrease  (cocl)  with  time.  The  two  morning 
to  evening  forecast  cases  did  not  reflect  sufficient  daytime 
heating,  even  with  clouds  present.  To  evaluate  the  pcor 


TABU  XT 

August-September  Lifting  Condensation  Levels 


E  ?.TE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCSTU 

29 

00 

781 

- 

- 

- 

— 

12 

196 

503 

— 

- 

- 

2C 

00 

520 

234 

227 

— 

— 

12 

206 

— 

0 

78 

— 

21 

00 

406 

— 

0 

297 

12 

364 

— 

— 

206 

CITE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

FCST8 

21 

12 

364 

105 

- 

— 

- 

1 

00 

467 

227 

613 

— 

— 

12 

M 

— 

580 

101 

- 

2 

00 

331 

- 

- 

150 

— 

12 

205 

— 

- 

- 

0 

3 

00 

577 

— 

— 

- 

0 

DATE 

GMT 

OBS 

FCST9 

FCST10 

FCST11 

FCST12 

3 

00 

577 

237 

— 

— 

— 

12 

284 

122 

325 

— 

— 

4 

00 

496 

- 

452 

424 

— 

12 

452 

— 

- 

341 

0 

5 

00 

525 

- 

• 

0 

DATE 

GMT 

OBS 

FCST1  3 

FCST14 

FCST15 

5 

00 

525 

374 

— 

12 

310 

263 

394 

- 

6 

00 

— 

— 

186 

194 

12 

- 

- 

- 

0 

performance  in  potential  temperature  forecasts,  the  subsi¬ 
dence  rates,  specific  humidity  and  diurnal  radiation  varia¬ 
tions  were  compared.  When  the  forecast  was  for  a  morning  to 
evening  period,  excluding  the  two  previously  mentioned  fore¬ 
casts,  continued  coding  was  predicted  when  heating  should 
have  been  predicted.  This  result  was  independent  of  subsi¬ 
dence  rates  or  presence  of  clouds.  Therefore,  for  morning 
to  evening  forecasts,  radiation  effects  were  not  sufficient 
to  introduce  the  needed  daytime  heating.  Prom  evening  to 
morning  the  major  problem  was  too  much  cooling.  This 
happened  in  five  of  the  forecasts  with  another  five  fore¬ 
casts  shewing  good  cooling  rates  but  bad  forecasts  for  the 
24-hcnr  point  since  the  previous  12-hour  forecast  was  bad. 
The  geed  coding  rates  occurred  when  the  subsidence  was  at 
its  largest  compared  to  the  excessive  lowering  forecasts. 
This  would  indicate  possible  problems  in  the  radiation 
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TABLE  XVI 


August 

-September  Potential 

Temperatures 

E  AT  1 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FC  ST** 

29 

00 

17.8 

- 

- 

- 

• 

12 

13.  5 

16.6 

- 

- 

- 

20 

00 

16.8 

15.5 

10.4 

— 

— 

12 

13.8 

— 

7.6 

13*7 

— 

21 

00 

15.8 

- 

— 

11.5 

9.8 

12 

14.2 

— 

— 

— 

7.6 

E  ATE 

GMT 

OBS 

FCST5 

FCST6 

FCST7 

FCST8 

21 

12 

14.2 

12.3 

- 

— 

— 

1 

00 

17.  € 

13.2 

12.1 

— 

— 

12 

M 

- 

10.3 

13.8 

- 

2 

00 

16.  1 

- 

— 

12.5 

— 

12 

12.8 

— 

- 

8.6 

3 

00 

15.  4 

— 

— 

- 

9.7 

DATE 

GMT 

OBS 

FCST9 

FCST10 

FCST1 1 

FCST12 

3 

00 

15.4 

9.9 

- 

— 

— 

12 

13.  1 

7.8 

12.2 

— 

— 

4 

00 

16.  C 

— 

10.8 

10. 1 

— 

12 

14.  1 

— 

- 

7.9 

12.2 

5 

00 

16.8 

— 

— 

- 

11.2 

DATE 

GMT 

OBS 

FCST1 3 

FCST14 

FCST15 

5 

00 

16.8 

11.6 

- 

— 

12 

13.8 

9.7 

16. 1 

— 

6 

00 

— 

— 

12.8 

11.2 

12 

- 

- 

- 

8.7 

•  • 
segment  which  would  he  cooling  the  potential  temperatures 

for  this  time  of  day. 

Figure  4.4  shows  the  observed  and  forecasted  inversions 
for  this  period.  Fey  factors  affecting  inversion  heights 
during  this  period  were  diurnal  heating  variations,  surface 
winds,  and  clouds.  The  temperatures  showed  a  definite 
diurnal  pattern  throughout  the  period  as  did  the  lifting 
condensation  level.  The  inversion  displayed  a  diurnal  trend 
only  part  of  the  time,  29/1200  GMT  through  31/0000  GMT  and 
3/0000  GMT  through  the  remainder  of  the  period. 

The  low  subsidence  values  contributed  greatly  to  the 
high  inversion  height  values  forecasted  during  this  period. 
Of  the  remaining  forecasts,  those  with  heights  forecast  to 
decrease  were  evening  to  morning  forecasts  with  subsidence 
occurring.  The  observed  heights  increased  (except  for  the 
2/OOOC  GMT  value  which  remained  constant)  due  to  strong 
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TABLE  XVII 


August- 

September 

Moisture  Value 

c 

DATE 

GMT 

OBS 

FCST1 

FCST2 

FCST3 

FCST4 

2S 

00 

8.5 

- 

• 

- 

- 

12 

8.7 

9.0 

- 

- 

- 

3C 

00 

9.1 

9.7 

6.9 

— 

— 

12 

8.  8 

- 

6.6 

9.3 

— 

31 

00 

9.0 

- 

— 

9.3 

6.4 

12 

8.  3 

— 

— 

— 

5.8 

DATE 

GMT 

OBS 

PCST5 

FCST6 

FCST7 

FCST8 

3  1 

12 

8.3 

8.4 

— 

— 

— 

1 

00 

9.9 

8.3 

6.3 

- 

— 

12 

M 

- 

5.7 

9.3 

- 

2 

00 

9.6 

— 

— 

8.3 

- 

12 

8.2 

- 

— 

— 

9.1 

3 

00 

8.0 

— 

— 

- 

8.9 

DATE 

GMT 

OBS 

FCST9 

FCST10 

FCST1 1 

FCST12 

3 

00 

8.0 

6.7 

— 

— 

— 

12 

8.  1 

6.1 

7.4 

— 

— 

4 

00 

8.7 

- 

6.3 

6.1 

— 

12 

7.9 

— 

— 

5.5 

9.0 

5 

00 

9.0 

— 

— 

— 

9.  1 

DATE 

GMT 

OBS 

FCST1  3 

FCST14 

FCST15 

5 

00 

9.0 

6.9 

— 

— 

12 

8.  3 

6.4 

9.3 

- 

6 

00 

— 

— 

8.3 

7.5 

12 

• 

- 

7.0 

winds  (weaker  and  more  westerly  at  2/0000  GMT)  which  lifted 
the  inversion  due  to  trixing. 

At  3 1/1200  GMT  an  observed  height  decrease  occurred  for 
the  morning  to  evening  forecast  time  but  the  forecasts 
showed  rises.  Tc  explain  this  requires  consideration  of  all 
three  factors  mentioned  above.  The  observed  height  values 
decreased  due  to  the  increase  in  cloud  coverage  and 
continued  onshore  flew  whereas  the  normal  situation  is 
cffshcre  flew  due  to  diurnal  changes.  The  radiation  compu¬ 
tations  alsc  appear  to  have  been  responsible  for  errors  in 
the  fciecasts  for  this  period.  Radiation  was  not  the  only 
problem  since  erroneous  specific  humidity  and  vertical 
velocities  contributed.  There  was  land  influence  on  the 
parameters  involved  in  height  predictions,  as  in  the  ether 
cases . 
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height  (meters) 


I.  S 0MH1H1  AMD  GENEBA1  C0H1EHTS 


It  is  easily  seen  by  the  large  root  mean  square  seer os 
and  the  large  differences  between  forecast  ?.n d  observed 
trends,  depicted  graphically,  that  the  model  did  net  yield  a 
good  prediction  cf  inversion  heights  at  Vandenbsrg  AFB.  The 
model  also  did  not  yield  very  good  forecasts  of  the  ether 
parameters. 

The  discussions  present  an  explanation  of  the  probable 
components  where  errors  occurred.  There  are  certain  prime 
factors  which  cause  the  greatest  problem.  The  largest  error 
producer  is  the  fact  that  Vanbenberg  AFB  is  a  land  base 
which  strongly  experiences  diurnal  variations  in  temperature 
and  winds.  The  primary  predictive  values  of  the  model 
(subsidence,  entrainment,  radiation  and  presence  of  clouds) 
are  also  influenced  by  the  land  effect.  Daytime  heating 
which  increases  temperature  and,  hence,  dissipates  clouds  is 
not  accctnted  for  in  the  MA EL  model. 

Topography  and  certain  wind  directions  increase  the 
turbulence  in  the  mixed  layer.  With  mountains  north  as  well 
as  south  of  Vandenberc  AFB,  this  condition  can  result  when 
the  offshore  flow  dees  not  set  up  during  the  night  and 
strong  northerly  or  southerly  flew  continues.  This  leads  to 
deepening  cf  the  mixed  layer  and  can  overcome  the  effects  of 
subsidence.  The  result  of  this  would  be  to  lift  instead  of 
lower  the  inversion. 

It  is  concluded  that  this  study  yielded  conclusive 
evidence  that  the  MAEI  model  is  not  appropriate  for  use  at  a 
coastal  station  because,  even  though  coastal,  land  influ¬ 
ences  play  a  significant  role.  The  results  from  the  study 
by  Gleason  (1982)  reveals  definite  contrasts  with  this  one. 
He  had  good  results  for  San  Nicolas  Island  and  was  able  to 
show  that  the  Q-methcd  of  calculating  vertical  velocities 
was  best. 
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V.  SOBHARY  A  HD  CONCLUSIONS 


Accurate  predictions  of  variations  in  the  height  of 
atmospheric  inversions  are  a  necessity  in  proper  usage  of 
modern  weaponry.  This  is  because  of  the  impact  changes  in 
inversicr  heights  have  on  the  ducting  of  weapons  guidance 
system  or  cf  tracking  system  signals.  Atmospheric  parame¬ 
ters  are  continually  changing  and  those  parameters  which 
govern  inversions  can  easily  change  their  entire  profile  in 
a  matter  cf  hours.  Forecasts  of  changes  in  inversion  prop¬ 
erties  must  start  with  forecasts  of  those  parameters, 
namely,  temperature  and  specific  humidity.  To  do  this  the 
marire  atmospheric  boundary  layer  model  used  at  the  Naval 
Postgraduate  School  was  utilized. 

Data  periods  which  covered  each  season  were  chosen  for 
the  Vandenbsrg  AFB  area.  All  four  cases  were  dominated  by  a 
high  pressure  system  over  the  eastern  North  Pacific  Ocean. 
Clear  as  well  as  clcudy  days  were  considered  to  allow  full 
range  testing  for  a  station  strongly  influenced  by  land 
conditions.  Data  for  these  periods  were  from  rawinsonde 
cbser vations,  satellite  information,  and  surface  and  500  mb 
synoptic  analyses.  Model  input  parameters  were  from  the 
rawinscrde  observations  and  satellite  derived  sea-surface 
temperatures. 

Results  from  each  model  run  were  compared  with  observed 
values  tc  determine  accuracies  of  the  model.  It  was  found 
that  the  model  did  rot  yield  accurate  predictions  during 
mcst  cf  the  cases.  Final  interpretations  of  these  inaccura¬ 
cies  always  point  to  strong  land  influence.  Diurnal  varia¬ 
tions  influenced  by  the  land  are  much  greater  than 
variations  over  water  and  are  not  reflected  in  the  model 
calculations. 
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He  end  result  cf  this  study  is  that  this  MAEL  modal 
dees  r.ct  *ork  for  Vancenberg  AFB .  It  is  not  currently  valid 
for  stations  strongly  influenced  by  land  induced  diurr.al 
variations.  The  re  develop  meat  cf  the  modal  to  tr-at  suer; 
factcxs  is  required  tc  properly  handle  the  land  influences. 
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